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A framework to calculate two-particle matrix elements for fully antisymmetrized three-cluster 
configurations is presented. The theory is developed for a scattering situation described in terms of 
the Algebraic Model. This means that the nuclear many-particle state and its asymptotic behaviour 
are expanded in terms of oscillator states of the intra-cluster coordinates. The Generating Function 
technique is used to optimize the calculation of matrix elements. In order to derive the dynamical 
equations, a multichannel version of the Algebraic Model is presented. 



I. INTRODUCTION 



Since 1980, the so-called Algebraic Model (AM) of the Resonating Group Method has been used in the investigation 
of bound and continuum states of nuclear systems. Initially the AM was applied to binary cluster configurations E-5|. 
Later on it was extended to describe binary clusters coupled to collective (quadrupole and monopole) channels |g-3|. 
Quite recently, three-cluster configurations were considered in the AM framework [p[~^l| . Such configurations play a 
significant role in light nuclei, in particular for reactions of astrophysical interest. 

The Algebraic Model represents the nuclear many-particle wave functions through their expansions in harmonic 
oscillator eigenstates. The use of a basis of square integrable states reduces the Schrodinger equation to a matrix 
equation. This procedure is well-known for bound states, but is also applicable to continuum states when the appro- 
priate boundary conditions are imposed on the expansion coefficients |I|P,|l2|-|l4[ . Thus the Algebraic Model provides 
a unified approach to bound and continuous spectra based on familiar matrix techniques. 

Solving an AM scattering problem requires two major steps. The first is to define the basis states relevant to the 
scattering channels being investigated and to compute the Hamiltonian matrix elements in that basis. The second 
is to consider the description of the asymptotic region and boundary conditions in the basis and to solve the matrix 
equation subject to those boundary conditions. 

We elaborate on both steps in this paper for a fully antisymmetrized three-cluster system, and introduce the proper 
three-cluster continuum boundary conditions. Such a description is still lacking in the literature. Three-cluster systems 
have been considered in some calculations, but there the continuum was approximated through sets of two-cluster 
configurations |l5^|l7|]. 

In our calculation the Pauli principle will be treated rigorously by taking full antisymmetrization into account 
for all nucleons in the system. The basis states will be defined by a composition of three s-clusters with "frozen" 
internal structure. The relative motions of the clusters are the relevant degrees of freedom. They are described by 
the Hyperspherical Harmonics method (HHM). It makes a natural formulation of the proper boundary conditions 
possible. The matrix elements of kinetic and potential energy in the three-cluster basis states are calculated using 
the Generating Function technique J|,[l8],[l9| , and methods for an explicit derivation will be presented. 

The asymptotic behavior of relative motion of clusters is obtained by considering large inter-cluster distances where 
the interaction between clusters approaches zero. We use the folding approximation to determine the asymptotic 
wave-functions. The folding approximation amounts to solving the scattering equations without nucleon-nucleon 
interaction and without inter-cluster antisymmetrization. 

In the case of neutral clusters, i.e. when inter-cluster Coulomb interactions are absent, the asymptotic equations 
are uncoupled in the hyperangular momentum K that is typical in the HHM. We will therefore introduce individual 
scattering channels characterized by quantum number K. In configurations with charged clusters the Coulomb 
interaction also couples if -channels in the asymptotic region. Again the folding approximation will be used to arrive 
at the appropriate asymptotic wave equations. A multiple channel approach which couples all i\-channels has to be 
used to solve the full AM system of equations. This extends the theoretical foundation of the AM given in [^0| for 
two-cluster systems. 
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In part II of this paper, the techniques discussed here will be applied to the three cluster configurations a + n + n 
for 6 He, and a + p + p for 6 Be. 

II. THREE- S-CLUSTER SYSTEMS 
A. Wave Functions 

The many-particle wave function for a three-cluster system of A nucleons (A = Ai+ A 2 + A3) can be written, using 
the anti-symmetrization operator A, as follows 

* (q l5 .., cu_i) = A [tf 1 {At.) * 2 (A 2 ) * 3 (A 3 ) V R (R)\ (1) 

where the centre of mass of the ^-nucleon system has been eliminated by the use of Jacobi coordinates . The cluster 
wave functions \&j (Ai) 

* i (A i ) = * i (q«,..,qg_ 1 ) (i = 1,2,3) (2) 

represent the internal structure of the j-th cluster, centered around its centre of mass R;. In our AM study these 
cluster functions are fixed and they are Slater determinants of harmonic oscillator (Os)-states, corresponding to the 
groundstate configuration of the cluster {Ai < 4 for all i). The (R) wave function 

*fl(i?)=*fl(qp ) ,q^ ) ) =*fl( qi ,q 2 ) (3) 

represents the relative motion of the three clusters with respect to one another, and qi and q2 represent Jacobi 
coordinates. In Fig. [l] we indicate an enumeration of possible Jacobi coordinates and their relation to the component 
clusters. 

FIG. 1. Two configurations of the Jacobi coordinates for the three-cluster system a + N + N. 

The state (^) is not limited to any particular type of orbital; on the contrary we will use a complete basis of harmonic 
oscillator states for the relative motion degrees of freedom. Thus the full ^-particle state cannot be expressed as a 
single Slater determinant of single particle orbitals. 

An important approximation is obtained by breaking the Pauli principle between the individual clusters, but 
retaining a proper quantum-mechanical description of the clusters, which is described by the wave function 

* F (qi,..,cu_l) = *i(Ai) ^ 2 (A 2 ) * 3 (^3) *r(R) (4) 

Because each cluster wave function is antisymmetric (they are Slater determinants) one is indeed neglecting the inter- 
cluster anti-symmetrization only. This results in what is known as the "Folding" model. It has the advantage of 
preserving the identities of the clusters and, if the intra-cluster structure is kept "frozen" , it reduces the many-particle 
problem to that of the relative motion of the clusters. 

The folding approximation will be the natural choice for calculating the asymptotic behavior of the three cluster- 
system, i.e. the disintegration of the system in the three non-interacting individual clusters. This amounts to the 
situation that all three clusters are a sufficient distance apart and intra-cluster interactions are no longer in force. 

The folding model is however also an acceptable approximation in the interaction region and can serve as comparison 
to the fully antisymmetrized calculations. In the folding model, the clusters interact through a local, inter-cluster 
potential called the folding potential. As in the current paper (Os)-determinants (Ai) are used to describe the 
internal state of the clusters, the folding potential will be easily calculated and is a sum of three terms 

VW = V (F) (Ri2) + V ( - F \R 23 ) + F (F) (R-3i) (5) 
where each term is simply the integral 

V^CRrv) = ]T I dr T dT v \y T (A T )\ 2 V(n -Tj +R™)| *^)| 2 (6) 

ieA T jeA v J 
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The coordinates R TlJ are associated with the relative position of the clusters 



1 

A, 



V 



(7) 



and sum to zero; they are equivalent to the qi,q2 Jacobi coordinates introduced earlier. In this way the folding 
approximation turns the three-cluster problem into an effective three-particle problem for the relative motion coordi- 
nates. 

Because the cluster states are fixed and built up of (Os)-orbitals, the problem of labeling the basis states with 
quantum numbers relates to the inter-cluster wave function only. This holds true whether one uses the full anti- 
symmetrization or the folding approximation. In a two-cluster case, the set of quantum numbers describing inter- 
cluster motion is unambiguously defined. In a three-cluster case, several schemes can be used to classify the inter- 
cluster wave function in the oscillator representation. In [p|jic[| three distinct but equivalent schemes were considered. 
One of these used the quantum numbers provided by the Hyperspherical Harmonics (HH) method (see for instance 
pl| §3[). This is the classification that we will adopt. Even within this particular scheme there are several ways to 
classify the basis states. We shall restrict ourselves to the so-called Zernike-Brinkman basis |24|. This corresponds to 
the following reduction of the unitary group U(6), the symmetry group of the three-particle oscillator Hamiltonian, 



This reduction provides the quantum numbers K, the hypermomentum, n, the hyperradial excitation, li, the angular 
momentum connected with the first Jacobi vector, I2, the angular momentum connected with the second Jacobi vector, 
and L and M the total angular momentum obtained from coupling the partial angular momenta l\, h- Collectively 
these quantum numbers will be denoted by v, i.e. v — {n, K 1 (7]^ 2 )£Af} in the remainder of the text. 
There are a number of relations and constraints on these quantum numbers: 

• the total angular momentum is the vector sum of the partial angular momenta liand I2, i.e. L = li+ I2 or 
IJ1-J2I <L< h + h. 

• by fixing the values of l\ and I2, we impose restrictions on the hypermomentum K = I1+I2, Z1+/2+2, Z1+Z2+4, . . . 
This condition implies that for certain values of hypermomentum K the sum of partial angular momenta li + 1% 
cannot exceed K. 

• the partial angular momenta l\ and I2 define the parity of the three-cluster state by the relation n = (— 1)' 1+ ' 2 . 

• for the "normal" parity states ir = (— 1) L the minimal value of hypermomentum is K m [ n = L, whereas K m i n = 
L + 1 for the so-called "abnormal" parity states n = (— 1) L+1 . 

• oscillator shells with TV quanta are characterized by the constraint N = 2n + K. 

Thus for a given hyperangular and rotational configuration the quantumnumber n ladders the oscillator shells of 
increasing oscillator energy. 



The Algebraic Model implements a method to solve the Schrodinger equation for quantum scattering systems, in 
particular for nuclear cluster systems. It is based on a matrix representation of the Schrodinger equation in terms 
of a square integrable basis, usually Harmonic Oscillator states, and boundary conditions in terms of the asymptotic 
behavior of the expansion coefficients of the wave function. In this paper we restrict ourselves to a presentation 
tailored to the treatment of three-cluster systems. 

In the case of three-cluster calculations, one needs to determine a proper approximation for the wave function (^). 
Consider an expansion of the relative wave function 



U(6) D 0(6) D 0(3) (8) 0(3) D 0(3) 



(8) 



III. THE ALGEBRAIC MODEL 



A. Asymptotic solutions in coordinate representation 




(9) 
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with v = {n,K, (lil2)LM} and a complete basis of six-dimensional oscillator states. It covers all possible types 

of relative motion between the three clusters. 

To obtain the asymptotic behavior of the three-cluster system, we consider the folding approximation. The assump- 
tion that antisymmetrization effects between clusters are absent in the asymtptotic region is a natural one. The relative 
motion problem of the three clusters in the absence of a potential can then be explicitly solved in the Hyperspherical 
Harmonics (HH) method (see for instance |2l|-|23]| ) . It involves the transformation of the Jacobi coordinates qi and 
q2 to the hyperradius p and a set of hyperangles Q. The inter-cluster wave function in coordinate representation is 
expanded in hyperspherical harmonics H'S (f2) where vq has been chosen as a shorthand for {l\l2)LM. 

In the absence of the Coulomb interaction this leads to a set of equations for the hyperradial asymptotic solutions, 
with the kinetic energy operator as reference Hamiltonian 



2 



2m 



d 2 5 d K (K + A) 



dp 2 p dp p 2 



E \ R KtUo (p) = (10) 



The solutions can be obtained analytically and are represented by a pair of Hankel functions for the ingoing and 
outgoing solutions: 



l K+2 1 

H$L (k P )/p 2 



where 



2mE 



One notices that these asymptotic solutions are independent of all quantum numbers v®, and are determined by the 
value of hypermomcntum K only. 

When charged clusters are considered the asymptotic reference Hamiltonian consists of the kinetic energy and the 
Coulomb interaction: 



2m 



dp 2 



5 d 
pdp 



\IC\\ 



P A 



-E \pL{p)\\ =0 



(12) 



The matrix \\K.\\ is diagonal with matrix elements K(K + 4), and ||^ e //||, the "effective charge", is off-diagonal in 
K and (hh)- Different if-channels are now coupled. A standard approximation for solving these equations is to 
decouple them by assuming that the off-diagonal matrix-elements of ||^ e //|| are sufficiently small: 



2m 



dp 2 



5 d 
pdp 



K(K + 4) 



Z, 



eff 



P 



E \ R K>U0 (p) - 



(13) 



The constants Z e f f depends on K and vq and all parameters of the many-body system under consideration. We will 
restrict ourselves to this decoupling approximation, but it is to be understood that its validity has to be checked for 
any specific three-cluster system. 

The asymptotic solutions then become 



R 



(±) 



(P) 



Wi Vtfl (2ikp) I pi 
W-iw (-2ikp) /pi 



where W is the Whittaker function, p = K + 2 and r\ is the well-known Sommerfeld parameter 

m Z eff 



n 



h 2 k 



(14) 



(15) 



As r/ is a function of K, l\ and I2 through the parameter Z e jf, the asymptotic solutions will now be dependent on 
K and v$. 
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B. Asymptotic solutions in oscillator representation 



The Algebraic Model relies on an expansion in terms of oscillator functions, and the asymptotic behavior of the 
corresponding expansion coefficients c„. It was conjectured (see for instance 20 1) that for very large values of the 
oscillator quantum number n the expansion coefficients for physically relevant wave-functions behave like 



(16) 



where p n — \Z4n + 2K + 6 corresponds to the classical turning point, b is the oscillator parameter, and ip is the 
hyperradial wave function. 

In the case of neutral clusters this leads after substitution of the hyperradial asymptotic solutions to the following 
expansion coefficients c< 



(±) 



(i) 



H 



(2) 
K+2 



(kbp n ) 
(kbp n ) 



(17) 



This result can be obtained in an alternative way |2q , |26fl by representing the Schrodinger equation, with the kinetic 
energy operator T as the Hamiltonian to describe the asymptotic situation, in a (hyperradial) oscillator representation 



E (». 



'A)) 



T-E 



m, (K, v ) ) c; 



JC,VQ _ 



= 



(18) 



This matrix equation is of a three-diagonal form because of the properties of T and the oscillator basis. Solving for 
the expansion coefficients c^' v ° leads to a three-term recurrence relation 



rjiK,V K,l/Q 

■ L n,n-l c n-l 



-iK,u 



E) c. 



K,v 



r K,v K,v a 
-n,n+l c n 



" — 
+1 — U 



where 



(19) 
(20) 



The asymptotic solutions (i.e. for high n) of this recurrence relation are then precisely given by (|T 
When the Coulomb interaction is present we again apply (|l6|) to obtain 



J±)K {2ikbp n ) I \fp~; i \ 

n ~ V 1 W-i^ {-2ikb Pn ) 1^ ] 



(21) 



In this case the oscillator representation of the Schrodinger equation is no longer of a tridiagonal form, and cannot 
be solved analytically for the asymptotic solutions to corroborate this result. 

It should be noted, that the above elaborations are valid for relatively small values of momentum k and sufficiently 
large values of discrete hyperradius p n , or when 



-IT < 1 

Pn 



(22) 



which defines the "asymptotic regime" ; it shows that for any value of k one can find values for n where the relation 
is satisfied. 

As we consider an asymptotic decoupling in the (K, vq) quantumnumbers, one will deal with asymptotic channels 
characterized by the (K, vq) values. So only in the internal (or interaction) region will states with different K and vq 
be coupled by the short-range nuclear potential and the Coulomb potential. The three-cluster system can therefore 
be described by a coupled-channels approach, where the individual channels are characterized by a single if-value, 
and we will henceforth refer to these channels as "A'-channels" . 



C. Multi-channel AM equations 



In the current many-channel description of the Algebraic Model for three-cluster systems, the channels will be 
characterized by the a specific value of the set of quantumnumbers K,i/q, whereas the relative motion of clusters 
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within the channel is connected to the oscillator index n. We will use K henceforth as a corporate index for individual 
channels, and assume it represents all K, vq quantum numbers. 

The Schrodinger equation can be cast in a matrix equation of the form 



EE( n >* 



K' m 



H - E 



m,K') cZ =0 



(23) 



We will now use a representation of the dynamical equations presented in 
formulation of the problem, the expansion coefficients are rewritten as 



As we will consider an 5-matrix 



(24) 



where, for the current channel K, the c^ K are the so-called residual coefficients, the dn K are the incoming and 
outgoing asymptotic coefficients (valid for all n). The matrix element Sk^k describes the coupling between the current 
channel K and the entrance channel JQ. 

As shown in p^ , p^ , |l^| , the Gn K satisfy the following system of equations for a given channel K 



m=0 



Hn-E 



m,K) c^ K =P { ±)K S n , 



(25) 



Hq being the asymptotic reference Hamiltonian, which equals the kinetic energy operator for uncharged clusters, 
and the kinetic energy operator plus Coulomb interaction for charged clusters. The right-hand side features I3^ K 
which is a regularization factor to account for the irregular behavior of the Cq . This factor allows one to solve 

( p5| ) for all values of n. The value of (3q K can be obtained for both reference Hamiltonians (i.e. with or without 
Coulomb). The set of equations (E^) for the asymptotic coefficients can then be solved numerically to different degrees 

of approximation depending on the requested precision. The Cn have the desired asymptotic behavior (cfr eqs |l7j 
and|2l]) 

Substitution of ( |24| ) in the equations (|2^) then leads to the following system of dynamical equations for the many- 
channel system: 



EE<S* 



K' m 



H-E 



m,K' 



E* 



KiK' 



K' 



a(+) K ' x x i V KK'(+) 
Pa n ,0 6 K'K + V n 



where the dynamical coefficients Vn K , defined in pot], are given by 



(26) 



±)KK' 



E(»>* 



V 



m,K>) c^ K ' 



(27) 



This system of equations should be solved for both the residual coefficients and the S'-matrix elements Sk'k- 

To obtain an appropriate approximation to the exact solution of (p6|), we consider an internal region corresponding 
to n < N and an asymptotic region with n > N . The choice of N is such that one can expect the residual expansion 

coefficients {c^ K } to be sufficiently small in the asymptotic region. Under these assumptions ( pif ) reduces to the 

following set of N + 1 equations (n — 0..N): 



EE 

K< m<N 



n,K 



H-E 



m,K')c^ K '-J2SK iK ' W )K ' S n , <W + C*' (+) 



K' 



a(~) K X X rrKKi(-) 
= -pfr d n fi 6 Ki K - V n ,y ' 



(28) 



The total number of equations for a given entrance channel Ki amounts to N c h(N + l), and solv- 
ing the set of equations by traditional numerical linear algebra leads to the N C ^N residual coefficients 

■fc^ ^; K = K m ia...Kxaax\ n — 0..N — l| and N c h S'-matrix elements {SkiK\ K — K m i n ..K max }. The set of equa- 
tions has to be solved for all N c h entrance channels labelled by Ki. 
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D. Numerical solution and convergence 



The numerical solution of the AM equations crucially depends on a proper choice of N, distinguishing the internal 
from the external region. The determining factor in this is the form of the potential energy matrix elements which, 
contrary to the short-range coordinate character of the potential, can be of a slowly descending nature in n. If this 
is the case, a sufficiently large value of N has to be chosen. 

In the case of three-cluster systems it is known from literature ]2£| that the potential asymptotically behaves as 1/ p 3 
in the hyperradius, with a corresponding effect on the matrix elements. It will be shown later on that the asymptotic 
form of the effective potential in the current case follows this behavior. It is well known p9|-|3l|| that potentials with 
an asymptotic tail 1/p 3 dramatically change the phase shift behavior in the low energy region and that special care 
should be taken to get convergent results. This cannot always be obtained by merely choosing a sufficiently large 
value of N. 

In reference [^0| a numerical strategy to account for long tails in the potential in the AM was developed. It 
dramatically improves the convergence of the results with significantly smaller values for N. We refer to [g(J f° r 
further details. 



IV. THE GENERATING FUNCTION METHOD 



A. General principle 



In this section, the general principles for calculating matrix elements in a three-cluster basis will be explained. The 
two main quantities of interest are: the overlap matrix, and the hamiltonian matrix. The former is of importance 
because of the proper normalization of the basis states. The latter is decomposed into the kinetic energy operator, 
the matrix elements of which are obtained mainly by group-theoretical considerations, the potential energy operator, 
which in our case will be chosen to be a semi-realistic two-body interaction based on a superposition of Gaussians, 
and the Coulomb contribution. 

In this work matrix elements for two-body Gaussian interactions will be derived. From these matrix elements of 
other functional forms of two-body interactions can be obtained using Gaussian transforms. This latter procedure 
will be followed to calculate the matrix elements of the Coulomb interaction. 

The basic principle of generating functions is well-known from mathematical physics. A generating function or 
generator state depends on a parameter, referred to as the generating coordinate, in such a way that an expansion 
with respect to that parameter yields basis states as expansion terms. A familiar example are the single-particle 
translated Gaussian wave functions 



(r |R) = exp 



— r 

2 



V2 R • r - -R 2 



(29) 



with the translation parameter R acting as generator coordinate. The choice of parametrization of the genera- 
tor coordinate influences the quantum numbers of the individual basis states that are generated. In a Cartesian 
parametrization R = (R x , R y , R z ) one generates the familiar Cartesian (j)n a .{Rx)4>n y {Ry)4'n z {Rz) oscillator states. 
With a radial parametrization R = RR (where the inverted hat stands for a unit vector) the expansion yields 



(r|R)= J2^niR 2n 1 

n,l,m 



Yi m (R)0„; m (r) 



(30) 



An underlying mathematical connection exists between such expansions, group representation theory and coherent 
state analysis |3^,[53| . In the present work we exploit the generating function principle to facilitate the computation 
of matrix elements. The matrix element of any operator between generating states is a function of the generating 
coordinates on the left and right 



X( R, R') 



(r|R) X 0(r|R 



(31) 



Expansion of this function will yield the matrix elements between the basis states. They can be identified in the 
expansion by the appropriate dependence on the generator coordinates 



X(R,R') = J2 E MdM l 'ri? 2n+ 'i?' (2n ' + °^(R)^ m ,(R)(^ m (r) 



(32) 



ill m nlm 
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Of course, one is not required to expand with respect to all parameters at once. Elimination of the angular dependence 
first, yields a partial generating function for the radial matrix elements: 



x(R,R') = ££y^(R)y rm ,(R')x (R, R') 

Im I'm' 

X lm;Vm ,(R,R') =J2Y, Af "i Af "'l' R2n+lR ' i2n ' +l ' ) (<W r ) |x| <t>n'l<m>{*)) (33) 

n n' 

Such partial generating functions will prove to be particularly useful when we apply the generating function method 
to the three-cluster problem with six generator coordinates (corresponding to six degrees of freedom) and an extensive 
set v = {n, K, (?iZ 2 )LM} of quantum numbers. 

The calculation of the matrix elements with the generating function method is a two-step process. The first step 
is the calculation of the generating function for the operator involved. Usually this is accomplished with analytical 
techniques. The second step is the expansion of the generating function w.r.t. the generator coordinates. Several 
approaches have been used in this respect. Explicit differentiation is one of them. Using recurrence relations for the 



expansion terms is another one |34|. In any case, the work involved here is straightforward but extremely tedious; 
both approaches are best implemented using algebraic manipulation software such as Mathematica or Maple. 

In this paper we introduce a representation of the generating functions in a manageable form to obtain explicit 
matrix elements and their connecting recurrence relations. 

B. Three-cluster generator state 

The customary generator state for the inter-cluster basis functions is given by (in what follows we shall use small 
q for the Jacobi vectors and capital Q for the corresponding generating coordinates) 



* (qi, q 2 |Qi, Q 2 ) = exp j-i (q? + q|) + y/2 (Qi ■ qi + Q 2 • q 2 ) - ^ (Q? + Q 



(34) 



The choice of parametrization is linked to the basis states one intends to generate. Associated with our choice of 
basis (Zernike-Brinkman p4|), we introduce hyperspherical coordinates. The hyperradius and hyperangles, both for 
spatial coordinates and for generating parameters, are defined by: 



qi+qf) ?i = pcos0, q 2 =psmt 



R= \jQi + Ql, Qi =RcosQ, Q 2 = i?sin9, (35) 
Using these, one expands the generating function ( ^) in hyperspherical harmonic functions: 

*(qi,q 2 |Qi,Q 2 ) =X>.,(p,0,qi,q2) K (R, 6, Qi, Q 2 ) (36) 

V 

where the full set of quantum numbers v (introduced previously) is involved in the summation. The oscillator basis 
functions are 

Mftfl.qi.qa) =Xn,K p K exp{-p 2 /2} L«+ 2 (p 2 ) H^ h)LM (9,qiA2) (37) 
and the generator coordinate functions are 

5„ {R, 6, Qi, Q 2 ) = M n ,K R K+2n H« lh)LM (6, Qx, Q 2 ) (38) 
Here H denotes the hyperspherical harmonic function 

i4 li2)iM (e,Q 1; Q 2 ) =A^ lZ2)LM $^ 2) (e){ Y;i (Qi) x Yz 2 (Q 2 )} LM 

^ li2) (e) =(cose) h (sine/ 2 P|lj4 +I (cos 26) (39) 
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From (p8[), one easily deduces the procedure for selecting basis functions with fixed quantum numbers v = 
{n, K, (Z1Z2) LM}. One has to differentiate the generating function (K + 2n)-times with respect to R and then 
to set R — 0. After that one has to integrate over 8 with the weight $^- 1 ' 2 ' > to project onto the hypermomentum 
K; one has to integrate over unit vectors Qi and Q 2 with weights Y/ im i(Qi) an d Y/ 2 m 2 (Q2) to project onto partial 
angular momenta. The order of these operations is not important and is a matter of convenience for each specific case. 
However these calculations, in particular those connected with integrating over hyperangle 0, are extremely extensive 
and cumbersome. For this reason, we introduce a new generating function appropriate for three-cluster calculations. 
We start from the function below which depends on seven generating coordinates namely: 

*(qi,q 2 |e,Qi,Q 2 ) (40) 

= (l + e)- 3 exp j-ii-^ ( q 2 +q |)+ J^_( Ql . qi + Q 2 . q2 )_i_l_(Q2 + Q |)| (41) 

It was used previously in a different context j7jj]i],|35|] to describe the coupling between monopole and two-cluster 
degrees of freedom. In those cases, the parameter Q generates basisfunctions of inter-cluster motion while parameter 
e generates collective monopole excitations of the ^4-nucleon system. Here, we will modify the function somewhat 
and use it only for the inter-cluster motion. We exploit the redundancy in the set of generating parameters (seven 
parameters vs. six degrees of freedom) and the fact that all expressions up to now are valid for complex generator 
coordinates also. We restrict the moduli of Qi and Q2 and set 

Qi = S,Q 2 = -iS (42) 



When the complex conjugate version of (|40|) is used, e.g. in the calculation of matrix elements, (|42|) is also complex 
conjugated. 

We now consider the new set of generator coordinates e,S, Q x , Q 2 , and substitute (|42| ) in (^0|), where the inverted 
hats on Qi and Q2 again indicate the angular components of both variables. This leads to 



* ( q |e,s,Q l5 Q 2 ) = ( 1 + e r 3ex P j-^^ (qi + q!) + Y^(Qrqi-*Q2-q 2 )j 

where the weights associated with each basis function are given by 

M^,Qi,Q 2 ) =Ml^ )LM H n ,K H) i2 e" S K {Y h {Qi) x Y, 2 (Q 2 )} LM (44) 

These are of a simpler structure and easier to use than ( |38[ ) because through (^) the dependence on the hyperangular 
coordinate has been eliminated. 

The full generating function for the matrix elements X u y of operator X now has the following general structure: 

jr(e,S,0i,0 a ;e , ,S / > $i,$a) 
= E (*" 0>, °> fc> fc) |x| ¥„, (p, 9, qi , q 2 )) $* (e, S, Qi, Q 2 ) $ v > (e',S', Qi,Q 2 ) 



X v y N ( K lh)LA W$ l ' 2)L ' M 'Nr hK Nn',K> H)' 2 {if ^ ^ S K S' K ' 



{Y1AQ1) x Y i2 (Q 2 )} iM {Y/^Qi) x Y^(Q 2 )} i/M/ (45) 

where again the shorthand notation for the quantum numbers v — {n,K, (/iZ 2 )LM) and v 1 = {n',K', {l^l^L 1 M') is 
used. As explained before, we will also consider partial generating functions which have been reduced with respect 
to a subset of generating coordinates and their corresponding quantum numbers. Most often, we will use a reduction 
with respect to the angular momentum dependence i.e. 

= 2-s X (hl2)LM-(l' 1 l' 2 )L'M'{^,S,e l \S') 

(hh) LM^l'^L'M' 

{^(Qi) x Yz 2 (Q 2 )} LA/ {Y^Qi) x Y^(Q 2 )} L , M , (46) 







with the partial generating function 

X(hl 2 )LM;{l' 1 l' 2 )L'M' ( e i S, 6 , 5') = 

£ ^n.a^^^M^.n.ai/^^M'A^^^^^^^^Vn.^n'.if' H)' 2 S K S' K ' (47) 

generating matrix elements for specified (Zi/ 2 )£M and (l[l' 2 )L'M' only. 

The asymmetry in the treatment of the different quantum numbers is motivated by the methodology in which 
the matrix elements will be used. Indeed, all quantities in this paper are calculated in the context of the Algebraic 
Model. As explained earlier the spatial asymptotic behavior is mapped onto the asymptotic behavior of the expansion 
coefficients in the oscillator basis. As for fixed K the n quantum numbers ladder through the oscillator shells, they 
will be needed for sufficiently high values in order to properly describe the asymptotic region. 



V. MATRIX ELEMENTS IN THE FOLDING APPROXIMATION 



In this section we derive the generating functions for the overlap and hamiltonian in the folding approximation, for 
several reasons. The folding approximation is indeed the natural representation for discussing the asymptotic behavior 
of the three cluster system, as the antisymmetrization between clusters vanishes at large inter-cluster distances. The 
calculation of generating functions in this approximation is also illuminating for the subsequent derivation of generating 
functions in a fully antisymmetrized setting, as the principles are identical, but the implementation is more complex. 
Finally, the folding approximation provides an interesting model to discuss the importance of antisymmetrization in 
the interaction region. 



A. Matrix elements for the overlap 



The overlap of two generating functions of the form ( |40|) is easily obtained, and can be written as 

2 



/(e,Qi,Q 2 ;e',Q 1 ,Q 2 ) =A- 3 ex P <ji^ 



i=l 



(48) 



where 



A = 1 — ee' 

After substitution of in fl4g|), one obtains 

f S S" 1 f S S" 1 

I (e, S, Qi, Q 2 ; e', S' , Q[, Q 2 ) = A" 3 exp | — Q x • Q'A cxp | — Q 2 • Q 2 j 



(49) 



It is interesting to note that the arguments of the exponential factors are diagonal in the generator coordinates S 
and S". To obtain an expansion in terms of angular momenta, the well known relation 



exp{a-b} = 4nJ2u(ab) Y;*„(a) Y* m (b) 



(50) 



can be applied, where i;(a;) is the Modified Spherical Bessel function of the first kind. Substitution of ( |50| ) in (|4S|), 
and applying traditional angular momentum coupling techniques leads to 

7(e,5,Q 1 ,Q 2 ;e',5',Qi,Q 2 ) 

= (4^) 2 £ i h (^pj i l2 (?pj Yr imi (Qi) Yf 2m2 (Q 2 ) Yr imi (Qi) YtmM) 

l 1 ,m 1 ,l2,m 2 



= (4?r) 2 ih(^)hJ^){YiAQi)><YiAQ2)} LM {YtM)ytM)} 



(hl 2 )LM 



LM 



(51) 
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The reduced generating function (cfr. (^T\j) then becomes 



l -(l 1 l 2 )LM;(hl 2 )LM 



(e,S;e',S') = (47ryA-*i h [ — )i h l — 



SS' 



SS' 



(52) 



This reduced generating function (|52|) is diagonal in the partial angular momenta l\ and 1%, and independent of total 
angular momentum L, thus valid for all angular momenta L compatible with l\ and 1%. 

The matrix elements with quantumnumbers K (also n) can now be obtained through a standard procedure, e.g. by 
differentiating (^2|) with respect to S and S' (e and e' for n). In particular, as (|5^) depends on the product of S and S' 
only, the overlap is diagonal in K . Likewise, the dependence on e and e' appears as a product in the factor A, so that 
the overlap is again diagonal in n. The fact that the generated overlap matrix elements are diagonal is a confirmation 
of the fact that these matrix elements were generated in an orthogonal basis. The calculation of the matrix is however 
not unimportant, as it provides a straightforward way to obtain the norm of the generator coordinate basisfunction 



B. Matrix elements for the kinetic energy 



In order to calculate the kinetic energy one can use the properties of the oscillator basis. The matrix elements of the 
kinetic energy of relative motion of the clusters are related to those of the oscillator potential by the virial theorem 



N 
N±2 



Nj = { N 
N) = 



Vo 
N±2 



V ) = ^E N 



Vo 



N 



where 



(53) 



(54) 



and En = hw [N + 3] = hu> [2n + K + 3] is the oscillator energy of relative motion. 

One easily obtains a representation of the oscillator potential in the manifold spanned by the generating function 



(q 1 +ql)*(qi,q 2 |e,S , ;Qi,Q 2 ) 
(l + c) 2 # + (l + £ )S-4 + 3(l + e 



de 



dS 



(55) 



so that the kinetic energy operator Tr can be represented by 



Tn = —hco 
2 



~ (1 ~ e)2 & + (1 ~ e)5 ZS +3(1 ~ £) 



(56) 



In a wider context, one can associate this operator with the generators TZ^ + \ TZ^ ^ and H\°> 



e 2 —+eS— + 3e 



de 



dS 



(57) 



of the 5*^(2, R) group classifying the space of relative motion of three-particle states, and whose irreducible represen- 
tations are labelled by hypermomcntum K so that 



Vr = \t\w 
Tn = \ hw 



n (o) _ n (+) _ n {- 



(58) 
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n, K, i/ 
n — 1 , K , Uq 



Using again the shorthand notation vq — (l\l2)LM, matrix elements are then readily found to be 

Tr 
Tr 



n, K, v Q J = -huj [2n + K + 3] 
n,K,u } = -^huy/(n +l)(n + K + 3) 
n,K,Vo^ = — — hajy / n(n + K + 2) 



(59) 



The full kinetic energy of the three cluster system must include the internal kinetic energy of the clusters: T 
Tr + T c i . As we consider frozen s-clusters only, this contribution is purely diagonal and equal to 



(60) 



C. Matrix elements for a Gaussian potential 



For a Gaussian two-body interaction with strength Vq and range a 



V(Ti,Tj) = Voexp 



(61) 



the folding potential (||) for (Os)-clusters can be calculated analytically. The result obtained from (^) is again of a 
Gaussian form but with a modified strength and interaction length now also depending on the oscillator parameter b 



^)(R™)=U z 3 / 2 exp{-^R^} 



(62) 



A T A V 
A T + A t 



The matrix element of a Gaussian potential between two generating functions of the form (40) is not necessarily 



diagonal in terms of the chosen Jacobi coordinates. In the folding approximation one can however easily find a set 
Jacobi coordinates in which the potential matrix element is diagonal. Both sets of Jacobi coordinates will then be 
related by an orthogonal transformation. 

In what follows, we distinguish the two types of Jacobi coordinates as follows: the original coordinates are denoted 
by q and were introduced in section II A, and shown in Fig. [j]; the diagonalizing coordinates will be denoted by x. 
An explicit instance of the latter are easily obtained. Indeed, for any choice of two clusters i and j (i ^ j), k being 
the third particle, one obtains a system of coordinates x, uniquely defined by index k, as follows, 



Xl 

x 2 



AjTj 



A* -\- A. a 



(63) 



At Aj A). \Ai + Aj) . 

CU= ^ = A^4-' C ^= A i+ A j+ A k (64) 

Each such Jacobi coordinate system k leads to a diagonal representation for the potential energy between clusters i 
and j. As both Jacobi systems q and x are related through an orthogonal transformation, we can invoke the Raynal- 
Revai theorem pa] . The latter states that any orthogonal transformation of Jacobi coordinates leads to an orthogonal 
transformation of the Hyperspherical Harmonics ( p9| ) in the wave function (pTl), preserving the hypermomentum 
quantum number K 

^ 2)LM (0 q ,qi,q 2 ) - £ O^IH ( k iX2)LM (^^^) (65) 

Ai A2 

The are known as the Raynal-Revai coefficients ]3r| . 
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This transformation can then be used to obtain the matrix elements of the potential in the original set of coordinates 
q through 



n,K,(l x h)LM V ri , K' , (l[l' 2 ) LM 



(q) 



= E E °xX°Xik (n,K,(X 1 X 2 )LM V n',K', (X[X' 2 )LM 



Ai A2 A' A' 



(x) 



(66) 
(67) 



where the matrix element on the right hand side, which is diagonal in the x representation, can be calculated in a 
straightforward way. 

To obtain the latter matrix element we consider a generating function of an identical structure as (^|) , in which we 
replace the q coordinates and Q generator coordinates by x and X. A generating function for the two-body matrix 
elements is then easily obtained as a product of two integrals over Xi and X2 leading to 



V (e,S, X 1; X 2 ; e',S', x' l5 x' 2 ) = V pre exp j ^X a ■ X^j exp j ^X 2 • X' 2 

V pre = V (AA)- 3/2 exp {-^ [(£'S) 2 + (£S') 2 ] } 
£ = 1 + e, £' = 1 + e' 

a = 1 - ee' = e + e - 

A = A + 7^' 

_ b 2 1 
7 a 2 ci a 



(68) 



By using (|50|), one again eliminates factorized terms and sums to obtain the following reduced generating function 
(cfr. (|47j)), in complete analogy to the procedure taken for (52): 



V, 



(hl 2 )LM;(hl 2 )LM 



(e, S; e', S') = V pre i h — — U 2 — — 



SS' 



A 



SS' 



(69) 



As was the case with the overlap, the matrix elements of the potential in the x coordinate system are diagonal with 
respect to the partial angular momenta l± and l 2 and do not depend directly upon the total angular momentum L ; 
this is a direct consequence of the characteristics of the operator. Expression ( |69| ) is again valid for all values of the 
total angular momenta L that are compatible with the given partial angular momenta l\ and l 2 , and thus represents 
a generating function for matrix elements with specific total angular momentum L. 

One obtains matrix elements with specific K and n quantumnumbers through the standard procedures (differenti- 
ation, recurrence relations, ...). So for example one obtains 



n, K, (hl 2 )LM 



V 



',K',(hl 2 )LM) 

I (x) 



d 

dS 



K 



dS> 



R" 



V(; 1 / 2 )LM;(;ii 2 )iAi" ( e ' &l £ 'i S) 



S=S'=0 
e=e'=0 



(70) 
(71) 



D. Matrix elements for the Coulomb potential 



The Coulomb interaction in the folding approximation between two clusters with Z r and Z v number of protons is 
given by (i) as 



V^ F) (R TV ) = Z T Z v e 2 Y, E [ dT T dT v \V T (A T )\ 



1 



ieA T j£A v ' 



\Vi - Tj + R T „| 



*.(^)l 5 



(72) 



A straightforward calculation of its matrix elements is impractical and very tedious. By however using the following 
Gauss-transform 



- = —= I dx exp{— r 2 x 2 } 



(73) 
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one rewrites ( |72| ) as 

v£ F) (™) = y V V / dT T dT v \^ T (A T )\ 2 / dxexp{- (r< - r 3 + R TV ) 2 x 2 }\ ^ V (A V )\' 



2Z-j-ZvC 2 1 



<i~ /' '^v{-^ri 2 z 



where 



= (l + *7 2 ) \ i = 2-M™, 7 = bx 



(74) 



(75) 



and its matrix elements can be obtained by integrating matrix elements (depending on z) of the Gaussian potential. 
Introducing the integration variable s 



tl 2 



1 + <7 2 ' 



transforms (\I4h to 



Vp(rv) 



2Z T Z v e 2 If 1 .! f R 2 TV \ 2Z T Z v e 2 f ,R 2 v . 

7=i 7=1 dss 2 ex P\-T7V s ( = 7=i erl (TTT) 

s/nb 2y/iJ Q I tb 2 J 0F6 tb 2 



(76) 



(77) 



This form shows that the integration can be reduced to a finite interval. It also shows that the Coulomb interaction 
between clusters does not behave as as could be expected. For very large value of R TV however ([7?]) properly 
reduces to Zt ^ c , due to the asymptotic form of the error function. 

E. Asymptotic behavior of the potential contributions 

As the folding model is used for defining the asymptotic channels, it is clear that the asymptotic behavior of the 
potential energy matrix elements in this model will be of vital importance for the rate of convergence of the AM 
solutions. The effective potential in terms of the hyperradius p is defined by integrating the folding potential over all 
hyperangles as formally indicated by 



W( P ) = J2Wrv(p) = W KM , h ( P ) 



K,h,h 



K,h,l 2 



(78) 



Its asymptotic behavior is then obtained for large values of p. 

In the diagonal representation with Jacobi coordinates x the calculation of W(p) between one pair of clusters r and 
v amounts to a straightforward integration over the hyperangles. 

The Gaussian interaction then exhibits the following asymptotic behavior 



Wrv{p) « (-1) 



K+3 *? f ru~a> 2/1+3 



K+h-h+l 

K-h- h 



K+h-h+l 

K-h- l 2 



which indeed shows a worst-case behavior (for h — 0) of the form 1/p 3 as predicted in section 



[II D 



The analogous evaluation for the Coulomb interaction leads to the following (exact) expression 



W TV (p) 



J K,l u h 
P 



(79) 



(80) 



with 
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7 e// - 1 7 7 p 2 rrr~ 



2 



3 

B(K - h - n - m + -,n + m + l x + 1) (81) 

and B stands for the beta-function. 

These results corroborate the fact that special care should be taken to get properly convergent results, even more 
critically when a Coulomb contribution between the clusters is present. 

VI. MATRIX ELEMENTS WITH FULL ANTISYMMETRIZATION 

When considering full antisymmetrization between all particles of the three-cluster wave function the normalization 
of the basis states becomes a non-trivial problem. The overlap, i.e. the matrix representation of the antisymmetrization 
operator A, needs to be explicitly calculated. 



As was carried out in the previous section one would normally start from the generating state (40) to obtain 
generating matrix elements in terms of e and e', facilitating the treatment of the hyperangular coordinates. We 
propose an alternative representation for the scaled generating state ( f40| ) which is more suited to our calculations, 
and hereto introduce the following integral transformation 

* (q|e, Q) = J dkcxp{-fc 2 } * (q; Q + V27k) (82) 

This allows to obtain ( f40| ) by scaling on the generator coordinate only. In other words, generating matrix elements can 
be obtained with the simpler generating state (E9) and scaled later on, reducing effectively the calculational burden. 



A. Matrix elements for the overlap 

Because the individual cluster states are Slater determinants one can use the familiar determinantal formulae to 
calculate the generating function. Starting then from the generating state (^9|) for the single-particle orbitals one 
obtains 

I (Qi, Q 2 ; Q'x, Q' 2 ) = J2 dH ex P I E ^ ' Q \ ( 83 ) 

The coefficients and B^ as well as the number of the terms depend on the specific type of three-cluster 

configuration, viz. the number of nucleons per cluster and their spin-isospin quantum numbers. 

In order to reduce (|83| ) with respect to the angular quantum numbers it is profitable to diagonalize the forms in 
the exponentials of (]83|). This can easily be achieved by diagonalizing its 2 by 2 coefficient matrix . This again 
amounts to making an orthogonal transformation from the original Jacobi coordinates q to new Jacobi coordinates 
x. This orthogonal transformation will induce a corresponding transformation of the basis functions which can be 
handled by the Raynal-Revai theorem as discussed in the previous section. The block-diagonal form of an exponential 
term can be written as 

expj^A, X 4 -X^| (84) 

We now introduce the scaling on e and e' by carrying out the transformation ( |S2"| ) on both generator coordinates X 
and X' for every (block-diagonal) term in (p3|), leading to 



I(e,X 1 ,X 2 ;e , ,X' 1 ,X£) = (A x A 2 y 3/2 exp 




X, • X: + £ ( e 'X 2 + eXf ) 



(85) 
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where 



A,- = 1 - A 2 ee' 



In order to further reduce the generating function in the partial angular momenta l\ and I2 we substitute ( }42| ) in ( |8q ) 
/ (e, S, Xi, X 2 ; e , S', X' l5 X 2 ) 
= (AiA 2 ) -3 ^ 2 exp |^=^ [SV + S' 2 e] } exp {^SS'Zi ■ X'J exp {^SS'X 2 ■ X' 2 | (86) 



By using (|50|) this produces in complete analogy to (|52| ) a generating function for all total angular momenta L 
compatible with l\ and Z 2 : 



1(h,l 2 )LM;(l 1 l 2 )LM ( e : "5) e > £") 



(87) 



If Ai = A 2 , the overlap depends on the factor 5*5" only and is therefore diagonal with respect to hypermomentum 
K, as was the case in the folding approximation. In the general case one will have Ai ^ A 2 , and (|87|) generates 
nondiagonal matrix elements in K. Hypermomentum is thus no longer a good quantum number for three-cluster 
systems, contrary to the folding approximation. 

The generating function for the overlap with K = 0, and consequently for angular momenta l\ = 1% = L = 0, 
is obtained immediately by putting S — S' — 0. This simplicity is an indication of the very suitable form of our 
generating function. 

The matrix elements I n '^°' ' v ° can now again be generated by a standard procedure, such as differentiation or 
recurrence relations. 



B. Matrix elements for the kinetic energy 



The matrix elements of the kinetic energy operator can be derived without use of a generating function. One of 
the effects of the antisymmetrization operator is to mix basis states within a fixed oscillator shell. The "diagonal" 
(i.e. within a shell) and "off-diagonal" (within neighboring shells) kinetic energy matrix elements are easily found to 
be connected to the matrix elements of the overlap by 



2n + K + 3 + - (A - 3) 



tK,Uq \K' ,Vq 



&2n 



2n+K,2n'+K' 



where the restriction on the quantum numbers to remain on the same oscillator shell (2n + K = 2n' + K') has been 
accounted for. 



C. Matrix elements for a Gaussian potential 



If one considers a Gaussian form for the nucleon-nucleon potential, and calculates a generating matrix element for 
the interaction using Slater determinants with individual orbitals of the form (p9[) one obtains terms of the form 



V (Qi, Q 2 ; Qi, Q' 2 ) =V (1- (f 2 exp { -( 



-, 2 



+ Bij q, ■ q; 



where 



2b 2 
2b 2 + a 2 



(89) 
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and Vq stands for any of the even (V13 and V31) and odd (Vn and V33) components of the TVA-interaction. Again 
b is the oscillator radius and a the range of the potential well. The first term in the exponent of (^) contains the 
factor £ and the vectors C1Q1 + C2Q2 and C[Q' 1 + G 2 Q 2 . The latter have a simple meaning: they define the distance 
between of the two clusters which respectively contain one of the nucleons in the interacting pair. 
In the diagonal representation (88) becomes 



V(X 1 ,X 2 ;X' 11 X' 2 ) =V cxp{-C 



(G,X 4 + G'^ 



+ X,.XJ) 



(90) 



i=i 



and the coefficients G,, C[ and Gi, G\ are trivially related by the orthogonal diagonalizing transformation. 

In order to eliminate crossterms in Xi • X 2 and • X 2 we introduce an additional transformation through the 
integral identity: 



1 f°° 

c?Zexp{-Z 2 + 2i(Z.a)} 

J — OO 



(91) 



leading to the following integral form for the block 

y(x 1 ,x 2; x' 1 ,x' 2 ) 



0^72 ex P jll X i X * ' X *j / dZ exp |-^ 2 + 2i(J2 {Gi Z • X 4 + G- Z • X-)| 
Using ( |S2|) on both X; and X^, one obtains after integrating over the corresponding and k^ : 

V (e; X X , X; e'; X[, X' a ) = (A 1 A 2 )~ 3/2 exp j£ A (e' A,X 2 + e A.Xf +2X r XQ | 

-A J°° dZexp j-AZ 2 +2iC^(& Z-Xi + % Z-X^j 



(92) 



(93) 



where again Aj = 1 — A 2 ee' and 



Gi + e'G^Xi , G\ + eGiXi 



A, 



A, 



A = |l + E ^ [eG. 2 + e'Gf + 2«'A 2 G 4 ^] | 



To obtain a reduced generating function for specific angular momentum quantum numbers, we again use expansion 
formulae of an exponential in terms of spherical harmonics. For exponential terms with a real scalar product we use 
(|50|), whereas for exponential terms with an imaginary scalar product we use |?7[ 



exp{i ab} = 4TrJ2i l h(ab) Y* m (a) Y* m (b) 



(94) 



Ini 



where j/(a6) is the well-known spherical Bessel function. 
Applying both (Fol) and (pi) to <M) leads to 



V ( e ,X 1) X;c' > Xi s X;) = (4^) 6 (A 1 A 2 )- 3 / 2 exp^ { -|- (e' X 2 + e Xf ) 

2—1 V 7 

— / dZexp{-AZ 2 }P klM ^ k ,, k[ ,^ JJj fe , (2^A 4 Z)j fc ,, (2^A 4 'Z) 



(95) 
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Before performing the integration over Z, we again make the substitution (M3) in (B5h to obtain 



V(e,S,X 1> X;e , ,S",Xi > X a 
(4 



(Ai A 2 )- 3/2 exp | ( e ' S 2 + e5 /2) | 

£ ^-ntf-^' ifei ^55') ifc2 (^') 
1 f 00 

/ dZexp{-AZ 2 }j fc , (26SZ)j fc ,, (2£S'Z)j fe , H2£ 2 SZ)j fe „ P kMl , K , K ,y> (96) 



where: 



PkiM,K,K.K>K =II( Y fc«( X *)- Y *'( X *)) (y M (X)-y m (Z)) (y*»(4)-Y*»(X<)) 

i=l 

= X] '^'{fc 1 ,fc 2 ,fe' 1 ,fe^fe;',fe^}:{(i 1 ,i 2 )L;(i' 1 J^)L':L"} 

(Ji,b)L;(Ji,i£)£';£" 

({Yj^Xx) x Y (2 (X 2 )} L • {{Y,i(Xi) x Y< (X 2 )} L , x Yl-(Z)} J (97) 

The intermediate (but redundant) index L" is connected to the integration variable Z of the integral transformation 
( |9T| ) . Because of the orthogonality between spherical harmonics and of the scalar character of the potential operator 
L and L' will be equal, and L" — after integration over Z: we therefore anticipate by simplifying ( |97|) to 

Pk 1 ,k 2 ,k' 1 .k! 2 ,k' 1 ' ,k' 2 / ] T^{k 1 ,k 2 ,k' 1 ,k' 2 ,k' 1 ' ,k' 2 '}-,{(l 1 ,l2)Li(l' 1 ,l' 2 )L:0}~7^ = 

jYi^Xi) x Y« a (X 2 )} £ x {Y i; (Xi) x Y< (X 2 )}J (98) 
The expansion coefficients are then easily shown to be 

r{fci,fc3,fci,fcJ,fci',fci'};{(Ji,ij)ii(Ji4)i;0} - ' | ; x ^ L J °fc0fe0 

H (-l) fcl (2fc, + 1) {2k> + 1) {2k1 + 1) | < ^ } ^° 0fe<0 C| 0fe(0 C^ fe ,, (99) 



Under the anticipative assumptions L = L' and i" = 0, the integration over the angles Z in (96) is now trivial. 
The remaining integration over Z is easily done after substituting the power expansions for the Bessel functions j/(x) 
and ii(x). The final result provides a tractable though very bulky result for Vn^LM-JhMLU ( e i &i e \ that is not 
reproduced here, as it carries no further additional information. 

Again, by using the standard procedures such as differentiation or recurrence relations one obtains the effective 
matrix elements. In particular the form of ( f70| ) remains valid, though the analytic differentiation preferably should 
be performed within an algebraic package such as Mathematica or Maple due to the bulkiness of the formulae. 



D. Matrix elements for the Coulomb potential 

The matrix elements of the Coulomb potential can now be most easily obtained from the Gaussian results. We 
consider again the Gaussian integral representation 

1 2 f°° 

= / cfeexp{-x 2 (r< - r,) 2 } (100) 

|r,--rj-| VttJo 

We use (p6|), replace a 2 by 1/x 2 , and introduce a new integration variable t — 7/(1 — 7) where 7 = x 2 b 2 , to obtain 
symbolically 
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Vc (e,S, X 1; X 2 ; e',S", X l5 X a 

= —7=l dx F (e,>S , ,X 1 ,X a ;e / ,S",!X 1 ,X a 
V 71 " Jo v 

= iT7=/ (e,5,^i.^;^,^,^) (101) 

Mutatis mutandis we apply the integral transformation directly to the generating matrix elements 

V(iij 2 )iM;(i;iyiM ( e 7 S; e , 5") 

V(l 1 h)LM;(l' 1 l' 2 )LM ( e ' e '' 

2 f 00 

= -p / ^(h^LM-.d'^LM i e i S; e', S') 

v 77 Jo 

2 f 1 di 



3 l "' / '(;i;2)iM;(/' 1 i^)LM (e,S;e',S') (102) 



6\/tF7 (l-t)tt3 

This however leads to a very intricate evaluation of the integral, and thus for the reduced generating matrix element 
of the Coulomb potential. 

A better procedure is to generate the quantum numbers K and K' first, by differentiation of 
V{hi 2 )LM;{l'l')LM ( e 7 S; e', S') on S and S' , then setting S = S' = 0. This leads to reduced generating functions 
for each K an d jf ' now only dependent on e and e', which are of a much simpler form and allow for an analytic 
integration in (|l02j). The further derivation on matrix elements for n and n' is then straightforward. 



VII. CONCLUSION 



In this paper we presented a framework for a microscopic three-cluster model within the Algebraic Model for scatter- 
ing. It was shown that it is possible to obtain matrix elements for fully antisymmetrized three-cluster configurations, 
as well as a proper description for the the three-cluster continuum in terms of a hyperspherical description. The 
corresponding AM equations in a multichannel description were also introduced. In the current work the individual 
clusters were limited to contain only s-orbitals, thus reducing the mass of clussters to that of a four-nucleon system. 
The latter restriction is however not a fundamental one, and was taken to restrict the analytical and calculational 
burden. 

In order to prove the validity and feasibility of the current model we will apply it to two specific three-cluster 
configurations believed to be of importance to astrophysical physiscs, a + n + n for 6 He, and a+p + p for 6 Be. These 
results appear in part II of this paper. 
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